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C*~> \ We explore the possibility that the dark energy is due to a potential of a scalar field 

and that the magnitude and the slope of this potential in our part of the universe 
£NJ ■ are largely determined by anthropic selection effects. We find that, in some models, 

. the most probable values of the slope are very small, implying that the dark energy 

Q , density stays constant to very high accuracy throughout cosmological evolution. In 

other models, however, the most probable values of the slope are such that the slow 
roll condition is only marginally satisfied, leading to a re-collapse of the local universe 
on a time-scale comparable to the lifetime of the sun. In the latter case, the effective 
■ equation of state varies appreciably with the redshift, leading to a number of testable 

CO , predictions. 

o \ 

1 I. INTRODUCTION 

o : 

,13 ■ It has long been suggested that both the old fine-tuning problem of the cosmological 

constant as well as the puzzle of the time coincidence may find a natural explanation through 
anthropic selection effects, in scenarios where the dark energy density pr> is a random variable 
' QQSBHBHBBHQHBHHIlS. This possibility can be easily realized 

>• , in the context of inflationary cosmology, where the local value of po may be determined 
by stochastic quantum processes. These processes may lead to rather different values of po 
in distant regions of the universe, separated by length-scales much larger than the present 
Hubble radius. 



A simple implementation of this idea is obtained [y, by assuming that the dark energy 
is due to a scalar field <fi (different from the inflaton field) with a very flat potential V(4>), 
which has a simple zero at = 0o with a nonvanishing slope s = |V A/ (0o)| : 

l/(0) = - S (0-0o) + O[(0-0 o ) 2 ], (1) 

where we have assumed for definiteness that V'((po) < 0. All that is required is that the 
slow-roll condition 

\V'\<H 2 M P , (2) 
is satisfied for values of the potential in the relatively narrow range 

\V\ < 10 3 M 2 P Hl (3) 



Here H is the present expansion rate and Mp is the reduced Planck mass, and we are adopt- 
ing the convention that any contributions to the vacuum energy (such as a true cosmological 
term) are included in the definition of V((f>). Larger values of |V| are uninteresting, since 
they would severely interfere with structure formation and with the emergence of suitable 
observers. During inflation, the value of the scalar field <fi is randomized by quantum fluctu- 
ations, and after inflation it stays almost frozen due to the flatness of the potential. Thus, 
the local value of the dark energy density po ~ V(0) will vary from place to place, but it 
will stay almost constant in time. In this situation, the probability for measuring a par- 
ticular value of pD is determined by a combination of inflationary dynamics and anthropic 
selection effects. As we shall see in the next Section, this approach to the cosmological 
constant problems shows remarkable agreement with observations, even with the crudest of 
assumptions. 

The purpose of the present paper is to extend this analysis to scenarios where the slope 
s of the potential is itself a random variable. Like pn, the measured value of the slope could 
be determined by a combination of inflationary dynamics and anthropic selection effects. 
A very large slope would cause a big crunch much before any observers can develop. If 
the distribution which is obtained after inflation favors large values of s, then a value of 
the slope which marginally satisfies (J2J) could be the most probable one to observe (Til fill. 
Til ]. Marginal slow-roll entails the consequence that the effective equation of state depends 



appreciably on redshift, po = w s (z)p D} through a function w s which contains a single 
parameter: the value of the slope s in our region of the universe. Thus, the equation of 
state (and its time evolution) may ultimately be determined by the condition that galaxy 
formation and the emergence of suitable observers is marginally allowed before the big crunch 
happens. Some observational signatures of models with a marginal slope have been discussed 



in 
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In Section II we review the case of variable pn at fixed s. In Section III we discuss two- 
field models of dark energy, where both p^ and s are random variables. Our conclusions are 
summarized in Section IV. 



II. VARIABLE PD 

A. Prior distribution 

As mentioned in the Introduction, a theory with variable p^ can be obtained from a 
scalar field with a very flat potential, as in Eq. (pj). During inflation, the field <fi undergoes a 
random walk of step size 5<p ~ H for each time interval St ~ H~ l , where H is the expansion 
rate during inflation. The steps are taken independently on each horizon volume, and this 
leads to spatial variation of <fi. The potential is very flat, and appreciable spatial variation of 
Pd after thermalization will only occur on scales much larger than the presently observable 
universe. 

In the limiting case when the potential is absolutely flat, the rate of expansion of the 
universe does not depend on the value of the field <fi. Then, because of the Brownian motion 
of the field <fi during eternal inflation, the field takes all possible values with equal probability. 
In other words, the volume distribution of the field at thermalization does not depend on 
the value of the field and takes the form 



dV* oc dcj). 



(4) 



When one takes into account that the potential is not entirely flat, the situation becomes 
more complicated. The probability distribution acquires some ^-dependence, which may 
be sensitive to a particular choice of the measure of p robability in an eternally inflating 
universe. This is a rather delicate issue, see e.g. [IJIi^l, but the final results may not 
be very sensitive to it because of the extreme flatness of the potentials suitable for the 
description of dark energy. It has been argued in ^1 that for a particular choice of the 
measure, and provided that certain generic conditions are satisfied, the volume distribution 
of the field <ft at thermalization preserves the simple form (j3J) in the narrow range of anthropic 
interest (J3J). We shall return to this issue in a bit more detail in the next Section, where 
the case with several dark energy fields is considered. As we shall see, additional subtleties 
arise in that context which require further discussion. For the rest of this Section, we shall 
assume that we are indeed in the situation where the flat distribution (JH) is valid. 

From the end of inflation until the present time, the field is heavily overdamped and 
remains almost frozen, giving a nearly time independent contribution to po- Thus, the 
"prior" distribution for the dark energy density is given by 

Consider, for illustration, the simplest linear potential 

V((f>) = a<\> . (6) 

If the slope of the potential is sufficiently small (as in most of the models of dark energy), 

a < 1(T 12 X> (7) 

then, according to Eq. (J2J), the field <fi practically does not change during the last 10 10 
years, its kinetic energy is very small. Therefore at the present stage of the evolution of the 
universe its total potential energy V(<j)) acts nearly like a cosmo logical constant. 
Similarly, one may consider a model 

V{4>) = \p 2 <P 2 + p A . (8) 

Here p\ is a true cosmological constant, and p?pa < 0, so that it is possible to have \pd\ 
very small even if the constant \p\\ is very large. Eq. (j2J leads to the condition 

H < 10- 120 M||p A |-^. (9) 

Note that the bounds on a and p do not correspond to a fine tuning, but just to a strong 
suppression. Possible mechanisms that could naturally account for such small values of 
parameters have been discussed in 

The potential vanishes at 4>q = —2pa/p 2 , and it can be easily checked that in the vicinity 
of this point we have 

\V'(<P)\ = s[l + 0( PD /p A )}, (10) 

where s is the slope at <ft — 0o and po ~ V(<p). Since the true cosmological constant is 
expected to be large in absolute value, we have po <^ Pa, and thus V'(<p) ~ const in the 
range of interest, so the potential should be well approximated by a linear function (flj). 
Substitution into © yields 

dV* oc dp D . (11) 



This means that all values of the dark energy density in the range (j3J) are equiprobable a 
priori. 

A linear potential as a simple model for dark energy was first considered in p|. Later 
it has been argued that this form of the potential is generic in the narrow anthropic range 



It is instructive to compare these models to the more traditional models of dark energy, 
with potentials of the type e _c< ^ or Q" 13 , where Q is the quintessence field jlsf. Generically, 
the potential of the quintessence field contains also a cosmo logical constant Vo, which, a 
priori, can be arbitrarily large and can have either sign. Thus, these models do not solve 
the cosmological constant problem. They also do not solve the coincidence problem, unless 
one fine-tunes the parameters of the potential. 

Since the quintessence potentials become asymptotically very flat, the scalar field Q also 
experiences quantum fluctuations during inflation. Therefore, in these simple models, one 
should expect that the typical value of the quintessence field becomes indefinitely large in the 
process of eternal inflation. As a result, these models become completely indistinguishable 
from the theory with a simple cosmological constant Vq. 

One can solve all of these problems, including the cosmological constant problem, by 
adding the interaction of the field Q with the curvature scalar, £-RQ 2 , and by multiplying 
the quintessence potential by n , where is a massless filed [TtJ- In the simplest case n = 1 
the potential of the quintessence models [28| in the regime when Q changes very slowly 
start looking very similar to our simple linear model ((T}. Because of the new term ~ 4>Vo, 
one can solve the cosmological constant problem in the same way as in the model with the 
linear potential jf| . The modified quintessence model will have some features distinguishing 
it from the model but overall it will be much more complicated. 



The distribution V*(pd) cannot be interpreted directly as the probability for measuring 
a particular value of po- If \pd\ is too large, so that it dominates prior to the galaxy 
formation epoch to ~ 10 10 yrs, then it will preclude the very existence of observers, and 
hence will never be measured. In order to implement this selection effect quantitatively, it 
seems reasonable to assume that we are typical observers in the ensemble of all observers in 
our thermalized region. The probability for measuring a particular value of pr> can thus be 
taken to be proportional to the number of civilizations in the universe which measure that 
value of pd, and we have 0, llfil ] 



Here, n(pr>, M)dM is the number density of galaxies of mass in the interval dM which will 
ever form in regions where the dark energy density takes the value /?£>, and N civ is the number 
of civilizations per galaxy. As a rough approximation we may assume that the integral is 
dominated by giant galaxies like the Milky Way, with mass M ~ Mq ~ 1O 12 M , and that 
N civ does not depend significantly on p D . For M > Mq, we may take N civ to be proportional 
to the number of stars in the galaxy, or to the mass of the galaxy, N C i V (M) oc M. Thus, the 
probability for measuring a particular value of pn is proportional to the fraction of matter 
f(M > MgiPd) which clusters in objects larger than Mq in regions with this value of pr>- 
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B. Full distribution 




(12) 



In the Press-Schechter approximation for determining the fraction of clustered matter in a 
ACDM model, and restricting attention to positive values of p D , one finds 

dV{p v ) cc f(M > M G , p D )dp D cx erfc [My 1 / 3 ] dy. (13) 

Here, we have introduced the variable y which is linearly related to the dark energy density 
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where a (M G ) is the present linearized density contrast on the galactic scale, Q m0 + Q D0 = 1, 
Qpo is the present value of in our local region, and the function F(x) is given by 



F(x) 



1 + x 



X 



1/2 



dw 



W 1 /6(l + w )3/2- 



(14) 
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FIG. 1: Comparison of anthropic predictions with observations. The curves represent the bound- 
aries of the 68% (solid) and 95% (dashed) confidence level regions predicted by the distribution 
(J13|). The cross represents the values inferred from WMAP observations, with 2a error bars. 



The distribution (|TH|) predicts that .33 < y < 6.0 with 68% confidence and .043 < y < 16 
at 95% confidence level. The observed value of y, given by 



F(n m /n 
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(15) 



is thus expected to fall into these intervals at lcr and 2a confidence levels respectively. The 
boundaries of the intervals define curves in the Qdq — a^ plane. These curves are shown in 
Fig. 1, where instead of using the present density contrast a (M G ) on the galactic scale, 
we use the more familiar quantity a 8 . For a given value of flno, o"s can be obtained from 
a (M G ) if the cosmological parameters such as the spectral index of density perturbations 
n, the dimensionless Hubble constant h and the baryon fraction Qf, are known. For these 
parameters we have taken the central values given by WMAP, n = .99, h = .72 and Qf, = .047 
[29j. Also shown in the same plot are the values inferred from WMAP for a 8 and Udo, which 
fall well within the anthropic predictions at the 95% confidence level. Given the simplicity of 



the assumptions which have been made, the agreement between predictions and observations 
seems quite remarkable. 

For negative values of pr>, the structure formation proceeds as usual until the time tp ~ 
(GIpdI) -1 / 2 , when the matter density p m becomes comparable to \pd\- At the moment when 
Pm = \pd\, the universe stops its expansion and starts recontraction. The effect of a negative 
Pd is in many ways similar to the effect of the slope, and to avoid duplication, we shall not 
discuss it here. Interested readers are referred to Refs. [3,0, where it is argued that the 
probability for p D < is less or comparable to that of pn > 0. In the following discussion 
we shall focus on the positive values of pp. 

III. VARIABLE p D AND SLOPE 

In Ref. jl6} two of the present authors considered several predictions of the anthropic 
approach to the cosmological constant problems, including the one we have just discussed 
in the previous Section. Two other predictions were that the equation of state should be 
that of a cosmological constant, w — — 1, and that the universe would re-collapse, but not 
before a trillion years. The latter predictions were based on the premise that, in a generic 
model, the slow-roll condition (J2J) is more likely to be satisfied by excess, by many orders of 
magnitude, rather than marginally. This seems clear in a model such as (JBJ), where the slow- 
roll parameter is fixed, and a strong suppression of p is required by the constraint (jHJ). If the 
suppression is due to some symmetry, then it is natural to expect that this symmetry will 
make the potential flat by excess rather than marginally. Indeed, the symmetry only knows 
about microphysics, and a marginal value \p\ ~ 10 _120 MJ,|pa| ^ 2 would itself represent a 
coincidence (or tuning) which requires separate explanation. 

However, a marginal value of the slow-roll parameter may be obtained more naturally 
in models where the slope is itself a random variable, by invoking again anthropic selection 
effects. It was pointed out in jTfJ that if the prior distribution favors large 5, then the most 
probable values of s could be the marginal ones. However, the authors of [la] argued (as we 
will show, incorrectly) that the prior distribution after inflation necessarily favors small s, 
and concluded that the equation of state w — — 1 should be expected in the general case. 
Dimopoulos and Thomas jl8| suggested, on the contrary, that the prior should generally 
favor large s, but offered no explicit model to justify this claim. Here, we shall examine this 
issue in detail in the context of a specific model. 

A. Prior distribution 

A variable slope is easily obtained by considering a model where we have several fields 
(p a instead of just one, so that on the hypersurface V(<fi, y?2, <Pn) = the gradient s(<p) = 
\W\ depends on <p a . 

During inflation, quantum fluctuations cause a random walk in field space which covers 
a distance |A<p a | ~ H(Ht) 1 ^ 2 in time t, where H is the inflationary expansion rate. If we 
start with some probability distribution dV*(<p a ,ti) at some initial time ti, the random walk 
causes "diffusion" of probability in field space, which tends to flatten the distribution as 
inflation proceeds. Hence, if we neglect the effect of the fields <p a on the expansion of the 
universe, we should expect that the volume distribution at the time of thermalization is 



given by 

dV*{(fa) oc JJcf</v (16) 

a 

In general, however, the potential of the dark energy fields V(<p a ) does contribute to the 
expansion, producing a nontrivial dependence of V* on (p a . Even though V is very small 
compared with the inflationary energy scale (at least in the anthropically interesting range), 
it causes a "differential expansion" which may accumulate during many e- foldings, biasing 
the distribution towards values of the field where V((p a ) is larger. This is simply because the 
volume of regions where the potential is larger grows faster. This leads to a field dependent 
distribution, 

dV*(<p a ) =F(<p a ) ]Jd<p a . (17) 

a 

In a region of size Aip in the field space, the characteristic quantum diffusion time is of 
the order r q ~ (Aip) 2 /H 3 , and the timescale on which the differential expansion becomes 
important is r de ~ (AH)- 1 ~ HM 2 /AV ~ HM 2 /\V v V\Ap. Diffusion will make the 
function F(ip a ) very smooth (or nearly constant) on scales Aip smaller than the smearing 
scale 

/ HM 2 \ 1/3 

Ayw, ~ H > ( 18 ) 

which is obtained by setting Td e ~ r q ■ On larger scales, F will generally have a nontrivial 
dependence on <p a . 

In the case we discussed in Section II, where there is a single dark energy field 0, one 
can argue that since the anthropic range for po is rather narrow, the corresponding range 
of 4> is also limited, and may easily be smaller than the smearing range (fTSj) . This has been 



used |l6l 1 171 ] in order to justify the use of the flat distribution (j3J) under certain generic 
conditions. However, this justification becomes less clear when we have several dark energy 
fields, since the anthropic range does not necessarily correspond to a small compact region 
in the field space, and F can vary significantly along the non-compact directions. 

Not much can be said about F in general, since its form depends on the overall shape 
of the inflaton and dark energy potential. (For a given potential, and with additional 
assumptions about the measure, F can in principle be calculated bv solving a suitable 
Fokker-Planck equation in the formalism of stochastic inflation (221 . 123L 0|)- To simplify 
our subsequent discussion, here we shall restrict ourselves to the case where differential 
expansion is negligible in the field range of interest. This is achieved for instance through a 
potential of the form 

U{il>,<p a ) = U{il>) + f(,/>)V(ip a ), (19) 

where ip is the inflaton and f a are the dark energy fields. The function f(ip) is normalized to 
unity in the thermalized phase, so that the potential V((p a ) becomes the dark energy density. 
If /(VO is sufficiently small in the range of ip corresponding to most of the inflationary phase, 
then the differential expansion can be neglected and the prior distribution will take the form 
©. 

The minimal number of fields required to account for variation of pu and s is n = 2. In 
this case, s = (Vj + V 2 2 ) l/2 and 

dV* oc dcpid(f2 oc I J\~ 1 dp£>ds, (20) 



where the Jacobian J is given by 



J = d( y > s ) = s - l e ab V, a V bc V, c (21) 

and V a = dV/dtpa- 

A simple example is given by 

V = a<f> + U{ X ), (22) 

where (p a = {<p, x) and a is a constant. The properties of the corresponding probability 
distribution can be easily understood by expressing it in terms of the variables pu and x- 
This gives 

dV* oc dpodx, (23) 

where we have used the Jacobian d(V, x)/d((f), x) = dV/dcj) = a = const. Thus, for a given 
interval of po, all values of x are equally probable. If the range of x is unbounded, then the 
character of the distribution is determined largely by the asymptotic behavior of U(x)- If 
U'(x) is a growing function of Xi then at large x the slope is s ~ and the probability 

is dominated by large values of the slope. For example, for U(x) oc x 2 , we have s oc x, and 
all values of s are equally probable (for sufficiently large s). Alternatively, if U'(x) ^ at 
|x| — > oo, then the value s = a is favored. 
For a potential of the form 

V = a<f>+- X n , (24) 

n 

where a, b and n are constants, the slope is given by 

s= ^ + b 2 x 2(n-l)y/2^ (25) 

and the Jacobian ()21j) is 

2n — 3 

Jocx 2n ~7socs-V-a 2 )5(^T). (26) 

Note that this expression depends only on s, but not on p D , and thus the distribution 
factorizes: 

dV* = V*{s)dsdp D , (27) 

with 

o o 2n — 3 

P*(s) oc s(s 2 - a 2 ) 2 ("-i). (28) 

The slope fl23j) satisfies s > a, so for consistency with © we have to require that a < 
HqM p , and the same logic as above leads us to expect this condition to be satisfied by excess, 
that is, a HqM p . Then, away from a small range of s near zero, the distribution (}28|) has 
a power-law form 

V»{s) oc s^s^. (29) 

Different behaviors of V*(s) can now be obtained with a suitable choice of the parameter 
n. For n — > oo, V*(s) oc s _1 , so all logarithmic intervals of s are equally probable. (This 
distribution is also obtained if the power law bx n is replaced by an exponential function 
e bx .) For n > 1, the distribution favors large values of s, and if n is chosen close to 1, 
the probability growth towards large s can be made arbitrarily fast. With n > 1, the 
prior distribution is non-integrable at large s, but in the next subsection we shall see that an 



effective cutoff is introduced by the galactic density factor n(p D , Mq), so the full distribution 
(H2J) is integrable. 

For readers who are concerned about the appearance of non-integrable distributions, even 
at an intermediate stage of the analysis, we note that the divergence comes from \ ~~ > 00 
and does not occur in models where the potential has the form (ffifj) only in a finite range 
of x- I n fact, Eq. ()24)1 with n = 2 can be thought of as an expansion of a more general 
potential (|2*2*|) in powers of \ near an extremum of U(x)- 

Finally, for n < 1, small values of s are favored, and again, by choosing n close to but 
below 1, the distribution can be made arbitrarily steep. In this case, the distribution ()28|) is 
non-integrable at s — a. This singularity is not smoothed out by the galactic density factor, 
and the prediction of models with n < 1 is that s = a with probability P = 1. Since a is 
expected to be very small, this prediction is observationally indistinguishable from s = 0. 
Once again, the divergence can be cut off if the power-law form of the potential P^)) applies 
only in a finite range of X- 



B. Galactic density 

The probability distribution for measuring given values of pr> and s is given by a straight- 
forward generalization of Eq. (j!2J) . 

dV(p D , s) oc J dMN civ (p D , s, M)n{p D , s, M)dV*(p D , s). (30) 

If the prior distribution dV* favors small values of s, then we expect the fields (p a to be deep 
in the slow roll regime. In this case, N civ and n are practically independent of s, and we 
recover the results of Section II. 

Suppose now that the prior favors large values of s. To simplify the discussion, we shall 
assume that the prior distribution has a factorized form (|2*7j) with a power-law distribution 
for s, 

V.{s)vlsP, (31) 

(3 > —1. As before, we shall assume that the integral in (p3H|) is dominated by giant galaxies 
of mass M > Mq. Assuming also that for such galaxies 

N civ (p D , s, M) = Mu civ (p D , s), 

where v civ is the number of civilizations per unit mass, we have 

dV(p D , s) = V(p D , s)dsdp D oc u civ (p D , s)f(p D , s,M > M G )s l3 dsdp D . (32) 

Here, / is the fraction of matter which clusters in galaxies of size bigger than Mq. This frac- 
tion depends on the relative magnitude of three characteristic times: the galaxy formation 
timescale to ~ 10 10 yrs, the onset of vacuum-like dark energy domination to ~ (Gpo)^ 1 ^ 2 , 
and the recollapse timescale t s due to the slope of the potential (we shall estimate t s shortly). 
If s is so small that t s is the largest of the three times, then the growth of density fluctuations 
effectively halts at tp, and the co moving density of galaxies can be estimated as in Section 
II. No matter how small, the slope eventually causes the field to roll down to negative values 
of the potential, ending in a big crunch. In the contracting phase, the density fluctuations 
start growing again, and one might think that any galaxies that failed to form at t < to 



would form then. However, "galaxies" that form at this epoch are likely to be very differ- 
ent from what we call "galaxies" now. At tp < t < t s , the dark energy density remains 
nearly constant, while matter density p m decreases exponentially with time, so at t ~ t s it 
is suppressed by an exponential factor. Moreover, in the course of the recollapse, the energy 
of the scalar field grows much faster than that of matter, and p m / Pd is further suppressed 
j^. Hence, the contribution of nonrelativistic matter (like CDM or baryons) to the mass 
of bound objects formed during the recollapse is utterly negligible. 

If tn is the largest of the three times, then the exponential suppression period is absent, 
but the rest of the above discussion still applies, and for t^$>t s the universe becomes scalar 
field dominated. Even at the onset of recollapse, t ~ t s , galaxies as we know them may not 
be formed. In our part of the universe, structure formation effectively stopped at t ~ tp, and 
the existing structures evolved more or less in isolation. This may account for the fact that 
discs of giant galaxies take their grand-design spiral form only relatively late, at z ~ 0.3. 
On the other hand, in a recollapsing universe the clustering hierarchy only speeds up at 
t > t s , and quiescent discs which may be necessary for the evolution of fragile creatures like 
ourselves may never be formed. 

This discussion suggests that for tjy < t s , the fraction of matter clustered in galaxies can 
be estimated as 

f( PD , s,M> M G ) ~ f(p D , M > M G ) (t D < t s ), (33) 

where f(pu,M > Mq) = /(ro,0, M > Mq) is the same which we used in the previous 
Section, which can be read off from (fT3"|). In the opposite case, tjj > t s , we expect that the 
density of habitable galaxies does not exceed the density of galactic-size halos that collapsed 
prior to t s . To estimate this density, we shall need an estimate of t s . 
The field equation for tp has the form 

p + 3(a/a)p> = s, (34) 

where a(t) is the scale factor. As long as the dark energy is subdominant, we have a(t) oc t 2 / 3 , 
and the solution of (JMj) is p = p^ + ^st 2 , where we have imposed the initial condition 0-^0 
at t — > 0. The dark energy density is then 

p D (t) = l -(p 2 - sip = p D <® - l -sH 2 . (35) 

Assuming first that t s < t D , we can disregard p§ , and the recollapse begins when the second 
term in (|3~5j) becomes comparable to the matter density p m ~ 1/Gt 2 , 

t s ~ (M p /s)^ 2 . (36) 

Alternatively, if the recollapse occurs after dark energy domination, then, for po > 0, 
a(t) oc exp(H D t) with H D ~ (Gpn) 1 ^ 2 . (Note that for p D < 0, the regime t s 3> to does not 
exist.) The solution of (|3*l)l is then p = ip(°> + (s/3H D )t, and the dark energy density is 

PD® = P { $ ~ (37) 

This equation applies as long as po remains nearly constant. Recollapse begins when the 
second term in (|3Tj) becomes comparable to the first, 

H D p D p% 2 

ts ~ — ~M^- (38) 



The boundary between the two regimes is 



t a ~ t D : s ~ Pd/M p . (39) 

Now, it follows from (p5fi|) that the matter density at £ s is p m (t s ) ~ sM p , assuming £ 8 < tp. 
This suggests that the fraction of matter in habitable galaxies in this regime is bounded by 

f(p D , s,M> M G ) < f(sM p , M > M G ) (t D > t s ). (40) 

In the estimates below, we shall use the value that saturates this inequality. 

We note finally that the dark energy density po in Eqs. (JH2|) should be understood as 
the value of pn immediately after inflation, that is, the quantity denoted by p§ in Eqs. (jSHJ), 



C. The number of civilizations 

To estimate the dependence of the number of civilizations per unit mass v civ on pz> and s, 
we now have to consider the role of two other characteristic timescales: £* ~ (2—3) x 10 10 yrs - 
the timescale on which most of the main sequence stars believed to be suitable for life explode 



as red giants (see [Hi, |23J] for more discussion and references) , and tj - the characteristic time 
needed for intelligent observers to evolve, tj is not likely to be much smaller than £*, since 
then it is not clear why it took so long for intelligence to develop on Earth. Carter |34| has 
argued that tj 3> £*, since the coincidence tj ~ £* is unlikely, considering that the evolution 
of life and the evolution of stars are governed by completely different processes. Note, 
however, that some seemingly unlikely coincidences may occur due to anthropic selection, 
tG ~ tr> being one example. Livio |3^ has suggested a simple model illustrating how tj ~ £* 
could arise. In any case, it seems reasonable to assume that 

h > U. (41) 

The time t* exceeds t G by only a factor ~ 3, but it will help to clarify the following discussion 
if we proceed as though t* ^> t G . This is justified in part by the fact that we will be comparing 
densities, which depend quadratically on time. 
Observers can exist only in the time interval 

t G < t < min{i s ,t*}, (42) 

and since according to (}4*Tj) this interval is shorter than tj, the number v civ is suppressed 
by a certain factor. Assuming that the origin of intelligent life is due to a single and very 
infrequent random event which has some constant probability to occur per unit time, we 
have 

v civ oc min{t s , t*} - t G . (43) 

In practice, many steps are necessary for the development of intelligent life, some of them 
occurring much more frequently than others. Assuming that, out of the total number of 
steps, there are k of them with typical frequencies smaller than l/[min{t s , t*} — t G ], then 
Eq. (pi3|) should be modified to 

Veto oc (min{t„ £*} - t G ) k . (44) 



Eq. (jUJ) assumes also that the steps which are needed to generate intelligence will produce 
the desired effect regardless of their time separation. This will not be the situation if there 
are relatively frequent catastrophes which occur at intervals shorter than (min{t s , — to), 
and which are serious enough to erase memory of any previously achieved steps. In this 
case, a linear expression such as (f4*3*j) is more appropriate. 

Our current knowledge about the number of steps k is very poor, and opinions differ quite 
vastly. Carter has argued that the effective upper bound to the total time biological 
evolution can proceed on earth, t b , is likely to be in the range t b — t e < t e /k, where t e ~ t* is 
the actual time intelligent life has taken to develop on earth. This formula reflects the fact 
that if many unlikely steps are necessary, then we are likely to have exhausted most of the 
available time before the emergence of intelligent life. If k is large, Carter's formula seems to 
indicate that some catastrophe is awaiting right around the corner which will erase life from 
earth within a timescale of order t b — t e <^t*. Carter rejected this possibility, and concluded 
that k > 2 was not very likely. However, Barrow and Tipler [3| have argued that there is no 
reason to reject a relatively imminent doom caused, for instance, by some instability in the 
evolution of the earth's atmosphere. This may render t b — t e much shorter than t*, and in 
this case there is no reason to expect that k should be small. Cosmic doom of the type we 
discuss in this paper is yet another way of obtaining t b — t e <C t*, since t b < t s and t s can in 
principle be smaller than i*. Qualitatively, however, our results will not depend strongly on 
k, and for the rest of the discussion we shall just take k = 1. 

The right-hand side of Eq. (}43|) takes different forms, depending on the relative magnitude 
of t s ,tD and t*. For t s > v civ is independent of s and pp,, 

v C iv ~ const (t s > t*). (45) 

For tc <C t s < t*, v civ oc t s , and using Eqs. (}36|) . (}38j) . we have 

v civ oc s~ 1/2 (t s < U,t D ), (46) 

Vdv oc Pn 2 s~ 2 (t D <t s <U). (47) 

The boundaries between the different regimes are given by Eq. (jHHJ) and by 

t a ~t* ( PD <l/Gtl): s~M p /tl, (48) 



t a ~U (p D >l/Gtl): s~ p%\M p U)- 1/2 . (49) 

The corresponding areas in the s — pD plane are sketched (not to scale) in Fig. 2. In our 
approximation, the factor Ud V vanishes for t s < to- (In a more realistic treatment, the 
density of galaxies would not strictly vanish for small values of t s . Galaxies would still be 
formed at high peaks of the density field, but their number density would be exponentially 
suppressed.) The boundary t s ~ t G is homotetic to the boundary t s ~ £* in the s — p D plane. 
For s > po/Mp it simply corresponds to the vertical line at s ~ Mp/t^, and for s < pp>/Mp 
it follows the curve 

s ~ pT{M P t G y l l\ (50) 
The behavior of v C i V as a function of pp> and s is illustrated in Fig. 3. 
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FIG. 2: Regions in the s — po plane, illustrating the different behaviours of v c { v with s. 



D. Full distribution 

We can now outline the general features of the full distribution (J32*|) . The effect of the 
galactic density factor n(pD, s, Mq) is, roughly, to cut off the distribution outside the square 
region po < l/Gt 2 G , s < M p /t G (These boundaries correspond to to ~ to and t s ~ to, 
respectively.) The fall-off, however, is rather mild, and n(pu, s, Mq) extends significantly 
outside the square. The cutoff is sharper in the s direction, due to the rapid decline of n civ 
with s. In region I V oc s@, in region II V oc s^ - ^ in region III V oc s' 3-2 , and V ~ in 
region IV. The distribution V for (3 = (n = 2 in Eq. (|24p) is illustrated in Fig. 4. The 
"logarithmic" distribution spoV^po, s) is illustrated in Figs. 5 and 6, for (3 = and (3 = 2 
respectively. 

For — 1/2 < /3 < 1, the distribution spoV^pn, s) is peaked along the line separating 
regions I and III (which corresponds to t s ~ t*). As we cross this "mountain range", moving 
towards larger values of s, the probability function drops in region III. The range terminates 
at point A where regions I-III meet. For — 1 < (3 < —1/2, the range continues beyond point 
A along the boundary between regions I and II. 

The probability distribution for can be obtained by integrating over s, 

V(Pd) = J V(p D ,s)ds. (51) 

The character of this distribution depends on the value of f3. 

For (3 ~ — 1 or smaller, the prior distribution is peaked on a very narrow strip near s = 
in the s — po plane. This strip is well within the flat plateau in v civ at low s (see Fig. 
3). In this case, integration over s just produces a constant factor independent of po, and 
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FIG. 3: Sketch of v c i v as a function of s and pp. For definiteness, we have used t* = 3tc- The 
slope s is in units of po/M pi whereas pd is in units of pc 

the posterior distribution coincides with the distribution V(pn) which is obtained for a true 
cosmological constant (s=0), given in (fTHjl . For (3 in the range — 1 < (3 < 1, but not too 
close to —1, higher values of s come into play. As shown in Fig. 3, the plateau in v civ is 
broader in the s direction at high pp, and this feature is inherited by the function V(s,pn) 
(shown in Fig. 4 for (3 = 0). Thus, the integration over s biases the distribution ()51|) towards 
larger p D , relative to the case of a true cosmological constant. This effect gets stronger as (3 
increases. 

Finally, for (3 3> 1, the distribution is pushed to the largest possible values of s. For 
(3 = 2, the distribution spdV(s, po) is plotted in Fig. 6. The probability is concentrated 
between the curves t s ~ t* and t s ~ to- If the value of j3 is further increased, the distribution 
gets packed more towards t s ~ tc, corresponding to (foT)|) . Larger j3 means that larger s is 
favoured by the prior distribution, and from (foT)|) . this means that the posterior distribution 
is peaked at even larger p^. The logarithmic distribution PdP^Pd) is shown in Fig. 7, for 
different values of /3, ranging from —1 to 2. The peak is indeed found to shift to larger values 
of pd as we increase (3. 

Fig. 8 shows the 1-cx and 2-a bounds on the variable y = Pd/pg as predicted by (Jo"Tj) . 
as a function of the parameter j3 in the prior distribution. The central value inferred from 
WMAP observations is yo ~ 0.1, which for (3 > —.5 lies outside the 2-a confidence level 
region. However, two things should be noted before jumping to conclusions. First of all, there 
is a large uncertainty in the measured value of y in our region of the universe. For instance, 
assuming the WMAP central values for the spectral index n pa .99, the dimensionless Hubble 
constant h pa .72 and the baryon fraction Qf, pa .047, and taking into account the 2 — a error 
bars for <r 8 and Qd (depicted as a cross in Fig. 1), we find from (fTB]) that the observed 
value of y lies in the range .04 < y < .48. Second, we should take into account that the 
predictions represented in Figs. 1 through 8 refer to the value of the dark energy density pr, 
at some very early time, when the scalar field is still frozen by the cosmic expansion. For 
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FIG. 4: The distribution V(s, pn) as a function of s and pp, for (3 = 0. As in Fig. 3, pr> is in units 
of pc- For this and the following plots, we have taken pq = (4/3)Mpt^ 2 , where to = t rec a r eJ 2 . 
Here a rec is the density contrast on the galactic scale at the time of recombination (to is then 
the time it takes for the linearized density contrast to become equal to 1 in the absence of a dark 
energy component). With this choice, the variable y which we used in section II is the same as 
Pd/pg- 

values of s > pc/Mp, the initial value of po may be larger than the value at the present 
time by a sizable fraction. For instance, if the slope s is such that the kinetic and potential 
energies of the dark energy field are approximately equal today (which may be considered 
a rather extreme case, although still marginally consistent with observations 0]) one finds 
that the dark energy density at very early times had to be larger by roughly a factor of 2. 
Hence, the initial value of Pd/pg m our region of the universe, may well have been anywhere 
in the range .04 < y < 1. The upper bound in this range is compatible with a value of (3 < 1 
(but not much higher) at the 2 — a level. 

Given the rough nature of our model and the uncertainties in the data, our conclusions 
must be regarded as very preliminary. Nevertheless, our results do suggest that low values of 
(3 < —0.5 should give a better agreement with the data. This leaves two possibilities: either 
(3 < —1, in which case the probability is sharply peaked at s ~ 0, and all predictions of 
Ref .[Trill for models with a fixed slope remain in force, or — 1 < f3 < 1, and then the recollapse 
is most likely to occur at t s ~ t*. Since is comparable to the present age of the universe, 
we can already expect to observe the signs of the slowdown of the cosmic expansion. 

IV. CONCLUSIONS 

The possibility that the smallness of the observed effective cosmo logical constant, as 
well as the puzzle of the time coincidence, may be attributed to anthropic selection effects 
is rather tantalizing. To implement this idea, one assumes that pr, is a random variable 
which takes different values in different parts of the universe (this could be due to stochastic 



10 



FIG. 5: The distribution sp£>V(s, pu) as a function of s and pp, for (3 = 0. Same conventions as 
in Fig. 4. 
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FIG. 6: The distribution spz>V(s, pr>) for (3 = 2. The rest of parameter values and conventions are 
as in Figs. 4 and 5. 

quantum processes which took place during inflation). Observers cannot live in regions 
where pc is too large, since galaxies cannot form there, so we should not be surprised at 
the smallness of the observed pp. Also, a simple analysis suggests that most observers will 
find themselves in regions which marginally allow the formation of suitable structures (e.g. 
galaxies of the type where observers are most likely to emerge). This would explain the time 
coincidence. 
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FIG. 7: The distribution pd'P(pd) as a function of pd/pg, f° r different values of /3. The solid line 
represents the standard case where pr> is a true cosmological constant. This is recovered in the case 
that P = —1, since in this case the prior distribution has a non-integrable peak at s = 0, where 
the slope of the potential vanishes. Also plotted are the cases (3 = — .9, —.5, 0, .5, 1, 2 (dashed lines, 
where higher (3 corresponds to longer dashes). The peak of the distribution shifts to higher /)p as 
/3 is increased. 

A pressing question regarding this scenario is whether it is possible at all to check its 
validity. In Ref. [HJ, two of the present authors ventured a few generic predictions of the 
anthropic approach to the cosmological constant problems. In particular, following up on 
the work of Martel, Shapiro and Weinberg it was argued that the variable y defined 
in Eq. (JT3J) should be y > .07 with 95% probability. Assuming a COBE normalized 
scale invariant spectrum of density perturbations, together with existing estimates for the 
baryon density, the anthropic argument suggested that the vacuum energy density parameter 
should be somewhat larger than = .7, or that the dimensionless Hubble rate should be 
somewhat smaller than h = .7. In Section II we have updated the comparison of predictions 
with observations by using the cosmological parameters as obtained from WMAP. Fig. 1 
shows confidence level plots for fl^ and cxg, corresponding to the 1 — a and 2 — a anthropic 
predictions, together with the values inferred from WMAP. 

The agreement with current data is rather encouraging, and it will be interesting to see 
how it evolves as the level of precision increases. Note that the confidence level regions 
in Fig. 1 are rather broad. This corresponds to a genuine large variance in the cosmic 
distribution of po- Hence, one may be led to the conclusion that future observations will 
not bring much excitement, since the overall picture will remain qualitatively the same even 
if the observational error bars shrink by a large factor. Nevertheless, we should recall that a 
number of assumptions went into Fig. 1. For instance, we assumed that the spectral index 
for scalar fluctuations, the Hubble constant and the baryon fraction are given by the WMAP 
central values. We have also used w ~ — 1 for the parameter in the dark energy equation of 
state. If the values of these parameters turn out to be different, the curves in Fig. 1 may 
shift significantly, putting some pressure on the anthropic explanation. 

Also, we must consider the fact that we are quite ignorant about the conditions which 
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FIG. 8: l-o" (solid) and 2-a (dashed) bounds on the value of the variable y = pd/pg, a s predicted 
by (|51|) . as a function of the parameter (3 in the prior distribution. Like in the previous figures, 
po is the dark energy density at very early times, and pc is the matter density at the time of 
galaxy formation. The central value inferred from WMAP observations is y ~ 0.1, assuming that 
the dark energy density has remained approximately constant through cosmological evolution up 
to the present time. 

are needed for the emergence of observers, an obvious drawback of the anthropic approach. 
However, if we can encode some of this ignorance in a few unknown physical parameters, 
we are in a position where we can predict something about the values of such parameters. 
By comparing the theoretical anthropic predictions with observations, we can find best fit 
estimates for the parameters, which may hopefully be confirmed some day by independent 
means. We have assumed throughout this paper that observers emerge predominantly in 
giant galaxies such as the Milky Way. This may be a reasonable assumption, but it is not 
an established fact by any means. Had we assumed that observers emerge predominantly 
in smaller galaxies, which form earlier on, the agreement with the data would become much 
worse. This reasoning was used in Ref. ^| to argue that the conditions for observers 
to emerge will be found predominantly in giant galaxies which complete their formation at 
redshift of order z ~ 1, but not much higher. This prediction seems hard to check at present, 
but hopefully much more will be known in the not so distant future about the properties of 
galaxies to confirm it or dispel it. 

In this paper we have considered models where the dark energy is due to the potential 
energy of several scalar fields. In the case of a single field, one assumes that the potential has 
a simple zero at = 0o with nonvanishing slope s (observers necessarily measure field values 
close to 0o, due to anthropic selection). If the slope is such that the slow-roll condition (j2J 
is satisfied by excess, then the equation of state will be indistinguishable from that of a true 
cosmological constant. But if (J2J) is satisfied only marginally, then there will be substantial 
evolution of the equation of state parameter w(z) with redshift. In models where the dark 
energy field has several components, both the dark energy density pr> an the slope s of 
the potential become random variables which take different values in distant regions of the 
universe (separated by distances larger than the present Hubble radius). The function w(z) 



(which is entirely determined by po and s) will therefore be different in each one of these 
regions. 

It was argued in Ref. that in the case of a single dark energy field, the slow roll 
condition (J2J was likely to be satisfied by excess, by many orders of magnitude, rather 
than marginally This leads to the predictions that the equation of state of dark energy is 
Pd = ~Pd with very high accuracy, and that the local universe will re-collapse, but not 
before another trillion years. It was also claimed that in generic two field models one should 
expect that small slopes would be favoured by the prior distribution, leading to the same 
predictions as in the case of a single field. However, the latter conclusion was based on an 
incorrect analysis of the prior distribution, which we have amended in the present paper. 

The prior distribution for the fields at the moment of thermalization can be obtained 
in principle from the inflationary dynamics. In the case of a single field, the anthropically 
allowed range (J3J) corresponds to a rather limited region in field space, and one can argue 
(under rather mild assumptions) that the prior distribution of the field will be almost flat 
within that range (l7l Efil]. On the other hand, in the case where the dark energy potential 
involves several fields, the range (JHJ) may correspond to a non-compact region in field space, 
and the prior distribution may be slowly varying in the non-compact directions. In this 
situation, the determination of the prior distribution from first principles is technically far 
more involved (and may require in general some further assumptions about the choice of 
the measure in an eternally inflating universe). Nevertheless, as argued in Section III, a 
flat distribution in field space can still be expected provided that the dark energy potential 
is sufficiently flat during inflation, so that its effect on the expansion can be completely 
neglected. In this situation, we have shown that there is a class of models where the prior 
distribution favors small slopes (in which case the conclusions of hold) but there is an 
equally broad class of models where large values of the slope are favoured a priori. 

The measured value of s is restricted by anthropic considerations, since if it is too large, 
the local region of the universe re-collapses before any observers have time to emerge. In 
Section III we have attempted to quantify this selection effect, and we have obtained pos- 
terior probability distributions for p^ and s. The problem of estimating the abundance of 
observers in regions with given values of pD and s has been split into two parts. In Section 
III.B we have discussed how the abundance of suitable galaxies is determined as a function 
of pd and s, and in Section III.C, we have analysed how the number of civilizations may 
depend on these parameters. There is of course much room for improvement in these esti- 
mates, but even at the rough level at which they stand, they do illustrate the fact that a 
posterior distribution which favors a marginal slope can easily be obtained in models where 
the prior favors a large slope. In this case we find that the universe is likely to turn around 
into contraction on a timescale comparable to the lifetime of the sun. This is a quite excit- 
ing prospect since it may lead to a potentially observable time-dependent equation of state 

HEBE!. 

We finally comment on the string theory motivated picture of a "discretuum" of flux 
compactifications with different values of po [13]. Recent work indicates that string theory 
does admit vacua with positive po (HI and that the corresponding spectrum of po may be 
rather dense [36|, suggesting the possibility of anthropic selection 0, 37, |38|. We note, 
however, that a dense spectrum of possible values for p D is only a necessary, but not a suf- 
ficient condition for explaining the value we actually observe. The probability distribution 
V(pd) depends on the prior distribution V*(po), and in order to obtain reasonable agree- 
ment with observations, the prior should not be too different from the flat distribution (j3J. 



However, nearby values of in the discretuum picture correspond to very different values 
of the fluxes. The parts of the universe with different values of pr> will have very different 
evolution histories, and one might expect that their probabilities will also be rather different. 
The arguments we gave in Sections I and IIIA for a flat prior distribution do not apply to 
this case. Calculation of probabilities in the discretuum remains an important problem for 
future research. 
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